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Intracellular transport along microtubules or actin filaments, powered by molecular motors such as kinesins, dyneins 
or myosins, has been recently modeled using one-dimensional driven lattice gases. We discuss some generalizations of 
these models, that include extended particles and defects. We investigate the feasibility of single molecule experiments 
aiming to measure the average motor density and to locate the position of traffic jams by mean of a tracer particle. 
Finally, we comment on preliminary single molecule experiments performed in living cells. 

I. INTRODUCTION 

Living cell is a highly organised structure that constantly needs to move its constituent parts from one place to 
another. It is therefore provided with a complex and accurate distribution system: the cytoskeleton, the network of 
biopolymers (actin, microtubules and intermediate filaments) that gives the cell its structural and mechanical features, 
functions as road system for the transport of organelles and vesicles; the motion of these objects is entrusted to motor 
proteins moving along these filaments. These motors are enzymes that convert the energy obtained by hydrolysis of 
an ATP (adenosin-triphosphate) molecule into a work (displacement of their cargo). Myosin V (on actin filament) 
and kinesin or dynein (on microtubules) are well known examples of these so called processive motors [l], Q • 

It has been observed that these motors can act cooperatively or interact with each other giving rise to collective 
phenomena. In particular, in some phase of the cell cycle, the motors can be expressed in high concentrations: it seems 
therefore natural to investigate analogies and differences with the traffic observed in a city. Experimental techniques 
as single molecule and fluorescence imaging have just started giving some hints on the complex behaviour of this 
system, but are still far from giving quantitative description of traffic situations. While waiting for experimental data 
some theoretical models have been developed to physically describe intracellular transport. 

II. DRIVEN LATTICE GASES: MODELS FOR INTRACELLULAR TRANSPORT 

A simple model to capture the behavior of many motors on a filaments needs to include three fundamental features: 
(i) the motors move in a step- like fashion on a one dimensional tracks and bind specifically to the monomer constituting 
the cytoskeletal filaments; (ii) cytoskeleton filaments present a specific polarization: the chemical properties of the 
track guarantee that the motion is always directed towards one of the two ends of the filament; (iii) the particles 
move according to a stochastic rule, i.e. they move upon a chemical reaction, the hydrolysis of ATP, which occurs 
randomly with a typical rate. 

The Total Asymmetric Simple Exclusion Process (TASEP) is a stochastic process first introduced to describe the 
motion of ribosomes on mRNA substrate [3| and encodes all these features. It rapidly became a paradigm of non- 
equilibrium statistical mechanics and one of the few example of exactly solvable systems 0, HI • I n this model each 
particle occupies a site on a one dimensional lattice and advances stochastically and in one direction. The most 
obvious observable is the average density profile of particles along the lattice. 

The system with open boundaries where particles enter the lattice with rate a at one end and leave with rate (3 at 
the other, shows a non trivial phase diagram where three distinct non-equilibrium steady states appear: a low density 
phase controlled by the left boundary, a high density phase controlled by the right boundary and a maximal current 
phase, independent of the boundaries. 

In a first attempt to construct a minimal model for molecular intracellular transport, one needs to add to the 
TASEP the fact that the tracks are embedded in the cytosol with a reservoir of motors in solution. This allows 
the motors to attach to or detach from the track. This led to the construction of the TASEP with Langmuir (i.e. 
attachment/detachment) kinetics or TASEP/LK model [1,0] depicted in Fig.QJu According to this model, in addition 
to the TASEP properties, particles enter (leave) the system with rate uja (^d) also in the bulk. All over the lattice 
they obey exclusion: two particles cannot occupy the same site. 
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FIG. 1: Schematic representation of the TASEP with Langmuir kinetics in the case of (a) monomers and (b) dimers. The 
allowed moves are: forward jump (with rate r = 1), entrance at the left boundary (with rate a), exit at the right boundary 
(with rate (3), attachment (with rate u>a), and detachment (with rate ujd) in the bulk. 



According to the rules described above, the rate equation in the average density at site i = 2 . . . N, (rii) can be 
written as: 

d t (rii) = (1 - rii)) - (m (1 - n i+ i)) + [u>a(1 - (rii)) - u D (rii)] . (1) 

The first two brackets describe the average current and the second one the on-off kinetics (source and sink terms) . At 
the boundaries this equation reads: 

d x (n x ) = a (1 - (m)) - (m (1 - na)> , (2) 
d x (n N ) = (njv-i (1 - "at)} - (un) • (3) 

Equation |l| shows a non-closed hierarchy in the correlation functions (i.e. (n,) depends on (niUj) and so on): this 
suggests the use of the mean field approximation (mrii+i) ~ (rii) (^t+i)- In the stationary state, the leading term of 
the continuum limit (N — > oo, (rii) ~ * p( x )) °f Eq. ([I]) reads: 

-d x [p(i-p)] + [n A (i-p)-n DP } = o, (4) 

supplemented by two boundary conditions: p(0) — a and p(l) = 1 — 0. When the solution of Eq. ((4| cannot 
be matched continuously with the left and right boundaries, the density profile displays a localized discontinuity (or 
shock) in the bulk (Figs.|2]>c). This translates into the emergence of mixed phases that were not present in the simple 
TASEP. For some sets of parameters the phase diagram can exhibit up to 7 kinds of coexistence (see Figj2^,). 

A key point in the study of TASEP/LK is the introduction of a mesoscopic limit where local adsorption-desorption 
rates u>a,d have been rescaled in the limit of large but finite systems ua,d = &a,d/N, such that the macroscopic rates 
are comparable to the injection-extraction rates at the boundaries [6|. Far from being only a simple mathematical 
trick that allows the competition of the directed motion with the on-off kinetics, this limit captures the fact that the 
motors explores a significant fraction of the track before detaching. This is precisely the limit of highly processive 
motors that biologically motivated our studies. Surprisingly only in the mesoscopic limit the density profiles show the 
shock. 



III. EXTENDING THE MODEL: DIMERS AND DEFECTS 



A. Dimers 



Many processive molecular motors (kinesins, dyneins and myosin V) are composed of two heads that bind specifically 
each to a subunit of the molecular track. A natural extension of the previous model towards a more realistic one 
consists in introducing non-pointlike particles in the system such as dimers (see Ref. || and Fig. [TJa). There are 
several challenging aspects in this problem: (i) the TASEP of particles of size £ (or t- TASEP) is known to have a 
non-trivial current-density relation j(p), different from the model with monomers (a, l9l. [Toj : (ii) even the simple on-off 
kinetics of dimers exhibits non-trivial dynamics, for example the stationary state is reached from an empty system 
through a double step relaxation process (III . [l2j . The coupling of an equilibrium process with two intrinsic relaxation 
regimes (on-off kinetics) to a genuine driven process (the I- TASEP) suggests interesting dynamical phenomena likely 
to result in new phases and regimes. 

In absence of exact solution, the main challenge has been to construct a refined mean field theory, based on 
probability theory, and to prove the consistency of the approximation for the two competing process (TASEP and the 
on-off kinetics). We have solved the mean field equation within the mesoscopic limit in the stationary state: 
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FIG. 2: (a) Phase diagram of TASEP/LK for uja = wd- One recognizes seven phases: in addition to the TASEP low density 
(LD), high density (HD) and maximal current (MC) phases; there are four more coexistence phases, namely the LD/HD, 
LD/MC, MC/HD and LD/MC/HD phases. The shaded region highlights the LD/HD coexistence where a localized domain 
wall appears, (b) Typical density profile in the LD/HD phase and (c) the corresponding current profile. At the matching point 
x w between the left (j a ) and right (jp) currents a domain wall develops and connects the left (p a ) and right (pp) density profile. 

The complete equation is formally similar to the one for pointlike particles (Eq(3|) : the current term (first bit in square 
brackets) must be balanced by the on-off term (second part). 

Exploiting the analytical properties of the solution of the mean field equation and the local continuity of the current, 
we constructed the global density profile. We have performed extensive stochastic simulations and the agreement with 
the mean field solution is excellent. As in the case of TASEP/LK of monomers a new phase coexistence region appears 
for some parameters. The main effect of extended nature of dimers on the phase behavior of the system is related to the 
breaking of the (particle-hole) symmetry of the model. This does have quantitative but not qualitative consequences 
on the density profile and on the phase diagram, which remains topologically unchanged. The origin of the robustness 
of the picture found for monomers can be traced back to the form of the stationary density profile which depends 
exclusively on the form of the current-density relation and of the isotherm of the on-off kinetics. In both the monomer 
and the dimer case the current-density relation is concave and presents a single maximum, while the isotherm is 
unique and constant: these two features are enough to determine the topology of the phase diagram. This robustness 
suggests that the TASEP dynamics washes out the interesting two-step relaxation dynamics that characterizes the 
on-off kinetics of dimers: the non-trivial outcome is that, in these systems, the diffusion (yet asymmetric) always 
dominates the large time-scale relaxation. 

B. Defect 

Another question that often arises in the study of these models concerns the role of some kind of randomness: the 
motion of the particles can be altered by structural defects of the track or by microtubule associated proteins (MAP). 
A good modeling requires to introduce some sort of disorder in the system: site-related or particle related, quenched 
or annealed. While disorder in TASEP has been treated with exact and approximated methods [HI, [13, [l5| ,the role of 
bottlenecks was never coupled to the TASEP/LK. As a preliminary study of the role of quenched disorder, it becomes 
particularly interesting to investigate the influence of an isolated defect (i.e. point-wise disorder) on the stationary 
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properties of the TASEP/LK. In Ref. [l6| we extensively studied this model where the defect has been characterized 
by a reduced hopping rate q < 1 (see Fig. [3]). 

As a consequence of the competition between the TASEP and LK dynamics, the effects of a single bottleneck in 
the TASEP/LK model are much more dramatic than in the simple TASEP [I?) (or the TASEP for extended objects 
[HI), where a localized defect was shown to merely shift some transition lines in the phase-diagram, but do not affect 
its topology. Here, new and mixed phases induced by the bottleneck have been obtained. 

As a key concept of our analysis, we have introduced the carrying capacity, which is defined as the maximal current 
that can flow through the bulk of the system. In contrast to the simple TASEP the spatial dependence of the current, 
caused by the Langmuir kinetics, makes the carrying capacity non-trivial: the defect depletes the current profile within 
a distance that we called screening length. This quantity increases with the strength of the defect and decreases with 
the attachment-detachment rates. The competition between the current imposed at the boundaries and the one 
limited by the defect determines the density profiles and the ensuing phase-diagram. When the boundary currents are 
dominant, the phase behavior of the defect-free system is recovered. Also, above some critical entrance and exit rates, 
the system transports the maximal current, independently of the boundaries. Between these two extreme situations, 
we have found several coexistence phases, where the density profile exhibits stable shocks and kinks. Indeed, above 
some specific parameter values the phase-diagram is characterized by bottleneck phases. Depending on the screening 
length imposed by the defect, which can cover the entire system or part of it, different phase-diagrams arise. The 
latter are characterized by four, six or even nine bottleneck phases, which have been quantitatively studied within the 
mean-field theory introduced in SecHTJ in fact for entrance and exit rates that exceed the critical value q/(l + q) (for 
which the bottleneck becomes relevant) the system can be split in two sublattices that can be treated as independent 
TASEP/LK with effective exit/entrance rates at the junction (Fig. [3^,). 

Our results were checked against numerical simulations, which brings further arguments in favor of the validity of 
mean- field approaches for studying the TASEP/LK-like models (Figs|3j3-c). The somewhat surprising quantitative 
validity of this approximate scheme can be traced back to the current-density relationship, which is correctly predicted 
by the mean-field theory. 

Eventually, we think that this study showed clearly that the presence of disorder in the TASEP/LK model, even 
in its simplest form, generally gives rise to quite rich and intriguing features and should motivate further studies of 
more 'realistic' and biophysically relevant situations, as in the presence of clusters of competing defects or quenched 
site-wise randomness (where the motors are slowed down at several points in the system). 
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FIG. 3: (a) Cartoon representing the TASEP/LK with a bottleneck: above a threshold density the system can be divided into 
two lattices where the results of TASEP/LK are valid, (b) Density and (c) current in a case where the bottleneck is relevant 
and competes with the boundaries leading to the formation of a double domain wall, theoretical results (dashed lines) are 
superimposed to numerical simulations (bold line). 
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IV. TOWARDS THE EXPERIMENTS 
A. Tracer dynamics 

All the models described so far would be rather useless without the possibility to measure quantitatively the density 
and to localize the shock. To this aim, simple bulk fluorescence imaging would be difficult to apply to quantitative 
experiments: single molecule analysis are far more promising. We proposed hence a method, based on single particle 
tracking, to measure the density of the system 19]. The idea is to use the simple TASEP (exact) results to reconstruct 
the density from the velocity or the diffusion constant of the tracer. It is known that the average velocity v = 
is related to the density p through the relation p = 1 — v. Exact results on a ring show that the same relation holds 

for the diffusion constant D = - — —. — 1 — p [2fJ ■ 

Since the rates scales with the size of the system, according to the mesoscopic limit, we suppose that the influence 
of attachment detachment on this relation to be negligible and the density to be locally continuous. We simulated the 
system TASEP/LK and measured the position x as a function of the time for several particles, in order to construct 
the probability density function P{x, t, 0, 0) to find a tracer particle at site x after a time t from its entrance in the 
system at site 0. From this function we can measure the velocity and the diffusion constant. As shown in Fig. 01 the 
shock is localized within a 10% precision through both the indicators. While the density is well reconstructed from 
the information on the velocity (Fig. 2^,), the hypothesis p = 1 — D works well only in the low density phase before 
the particle arrives at the shock. Once the particle has passed through the shock the hypothesis on the continuity of 
the density breaks down and the relation is not valid anymore while the loss of particles in the high density phase 
makes the system subdiffusive. Yet, the quantity 1 — D can be used to localized the shock (Fig. |4Jd). 

Numerical simulations show that, on system of realistic size (i.e. lpm, roughly 100 sites), analysing the trajectories 
of a hundred particles is enough to reconstruct the density profile and localize the shock (after a time moving average 
resulting on a smoothing of the data). This suggests that single molecule experiments observing local features could 
give information on the global scale. 




FIG. 4: (a) Density profile of TASEP/LK derived from the analysis of 100 tracer particles: wiggly line match rather well 
the theoretical prediction (dashed line), (b) Superposition of the first two moments of P(x,t), (x) (t) (dashed) and V^ar[a;(t)] 
(wiggly) in time: according to the simple TASEP results the two quantities should be the same; this is true only until the 
particles have reached the shock, afterwards the variance deviates. 



B. Single molecule in the cell 

So far much information concerning processive molecular motors was provided by single molecule experiments. 
Typically a latex bead is attached to the motor to make it visible and to manipulate it: this allowed us to measure 
not only velocity and processivity of the motors, but even the forces they exert, the steps and (in some lucky cases) 
substeps [2lll22. 23, 24]. In- vitro observations can be quite controversial: some experiments did not show any difference 
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upon increasing the concentrations of motors, as if traffic effects were not relevant 25[. However the situation in the 
overcrowded environment of the cell is hardly reproducible in vitro and could reserve more surprises. 

Since single molecule experiments have often been criticized because of their distance from the real systems, many 
groups are cautiously moving the single molecule experiments directly into the cell [2(| [U [H, [2t| . In this situation 
the main problem is to have a good signal- noise ration which is hardly achievable with the usual fluorescent probes. 
Recently a method to mark and observe single molecule in vivo by using quantum dots (QDs) has been proposed 
[28l [30j . In contrast with the usual fluorescent probes, the QDs do not bleach, are excitable on all the visible 
spectrum and show a very narrow emission band. The only drawback is that they blink without a typical timescale: 
this inconvenient can be avoid by taking longer series of images. In Ref.[3(| kinesins are biotinilated to bind to a 
streptavidinated QD. The conjugated particles so obtained are introduced in Hela cell by pinocytosis followed by 
osmotic shock and the motion of the QDs is observed with a customized fluorescence microscopy setup and a fast 
CCD camera. The QDs signal on 25 pixels is fit with a Gaussian to obtain sub-pixel resolution (FIONA [3l|). The 
blinking of moving QDs is taken as an evidence that we are working in single molecule. The results on kinesin speed 
(570 ± 20nm/s) and processivity (1.73 ± 0.06s, i.e.^ 1/im in space) are compatible to the ones known from in-vitro 
experiments (Fig. EK-b) (34[ . 
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FIG. 5: Histograms of the speed (A) and processivity (B) of kinesin in Hela cells. The results confirmed in vitro measurements. 
(C) Picture of the trajectories (arrows) of a QD labeled Kinesins: image obtained choosing the maximal intensity over many 
frames of a movie with long integration time (100ms). 



V. CONCLUSION 

In this proceeding we have reviewed a model introduced to describe intracellular transport @ and its possible 
extensions towards a more realistic picture: we investigated how the introduction of dimers [1] and the presence of 
disorder on the track [l6| would affect the known results. These theoretical approaches while aiming to a realistic 
description of the traffic phenomena, are inspiring statistical models interesting in its own right (see e.g. [32I [33|). 
The latest experimental successes in tracking individual motors in living cells combined with an appropriate analysis, 
presented in the last section, could not only confirm the theories but also provide ispiration for a complete picture of 
the cell logistic. 
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